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We study the Kohn-Sham scheme for the calculation of the steady state linear response 
AnL 1 (r) cos u)t to a harmonic perturbation Ai/ 1 -^! - ) coswi that is turned on adiabatically. Al- 
though in general the exact exchange-correlation potential v xc (r,t) cannot be expressed as the 
functional derivative of a universal functional due to the so-called causality paradox, we show that 
for a harmonic perturbation the exchange-correlation part of the first-order Kohn-Sham potential 
cosujt is given by v2c (r) = 5K^c /SnQ (r). K^J is the exchange-correlation part of the 
second-order quasienergy K v . The Frenkel variation principle implies a stationary principle for 
Kv [tow ]- We also find an analogous stationary principle and KS scheme in the time dependent 
extension of one-matrix functional theory, in which the basic variable is the one-matrix (one-body 
reduced density matrix). 

PACS numbers: 31.15.ee,32.10.Dk,71.15.Mb 



I. INTRODUCTION 

The description of a time dependent quantum state is 
fundamentally different from the description of a ground 
state because there is not a minimum principle for the 
former. As density functional theory 1 (DFT) is based on 
a minimum principle, it was not obvious that it could 
be extended to time dependent situations. Yet such an 
extension, a time dependent DFT (TD DFT), was estab- 
lished by the Runge-Gross (RG) theorem^ which asserts 
that the time dependent density of a many-electron sys- 
tem determines the time dependent external potential up 
to an additive purely time dependent function, assuming 
a fixed initial state \&(io) — ^t - The theorem applies 
to both interacting and noninteracting systems. This is 
of great importance for applications because it implies 
that the density n(r,t) of an interacting system can be 
reproduced by a noninteracting system with an effective 
potential v s (r,t), provided the initial state $(to) = $t 
of the noninteracting system is chosen to be compatible 
with n(r,£).— The potential v s (r, t), which is called the 
time dependent Kohn-Sham 4 (KS) potential, is a func- 
tional of n(r, t), ^>t and $t . Its exchange-correlation 
part v xc (r,t) can be defined from the equation 

v s (r, t) = v(r, t) + v H (r, t) + v xc (r, t), (1) 

where vh(y,£) is the time dependent Hartree potential. 
In contrast to v xc (r) in static DFT, the exact v xc (r,t) 
cannot be expressed as the functional derivative of a uni- 
versal functional of n(r,t). This is a consequence of the 
so-called causality paradoxi^i^ In TD DFT, neither a 
general minimum principle nor even a stationary princi- 
ple has been founds Such a principle might be helpful in 
the search for accurate approximations to v xc (r,t). The 
quantum mechanical action principle does not lead to a 
stationary principle of the form 5 A = 0, where A is a 
functional of n(r,t), because its density-functional for- 
mulation contains boundary terms^ 



An important application of TD DFT is the calculation 
of dynamic polarizabilities and excitation energies i 9 ' 10 : 11 
These can be obtained from the frequency dependent lin- 
ear response function xi 10 ) = x( r > r 'i <^)- 32 In the time do- 
main, the retarded linear response function is defined as 
X(r, t, r', t') — 5n(r, t)/5v(r', t'). The KS system repro- 
duces self-consistently the linear response of the interact- 
ing system to a perturbation Sv(r, t). Thus, x( r ) t, r', t') 
is related to the KS response function Xs( r ,^ r ',£') = 
5n(r, t)/5v s (r', t 1 ) by the Dyson-like equation 1 ^ 

X(1,1')=X.(1,1') + / d2 d3 X .(l,2) 

x ( Wc (r 2 , v 3 )S(h - *s) + fxc(2, 3)) X (3, 1'), (2) 

where i = {r^U), ^(r^ry) = |r< - r^ -1 and f xc (i,j) = 
8v xc (i) 1 5n(j) is the exchange-correlation kernel. If the 
system is in its ground state at t = to, then x(^') de- 
pends only on the time difference t — t' and not t and t' 
individually. Thus, from the Fourier transform of ^ one 
obtains 

X(r,r';cj) = x s (r,r';w) + J d 3 r 2 d 3 r 3 x s (r,r 2 ;w) 

x (u c (r 2 , r 3 ) + /xc(r2,r 3 ; w))x(r 3 , r'; u) (3) 

if fxcifjj) exists. Gross and Kohn state that "For the most 
general situation, we do not know whether f xc exists..."— 
Eq. ([3]) implies the following formal representation: 

f xc (r,r';w) = xj l (r,r';w) - x _1 (r, r'; - u c (r,r'). (4) 

Therefore, f xc {u) exists whenever x( w ) an d Xs(w) are 
invertible. Although the RG theorem guarantees that 
X(t, t') is invertible (subject to the condition that the 
perturbation is analytic at t = to), the frequency depen- 
dent response function x{ u ) IS not always invertible.— 
As x( w ) is the Fourier transform of x(t — t'), one might 
ask why the invertibility of x{t — t') does not imply the 
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invertibility of x{u). We shall address this question in a 
later section. The problem of the invertibility of the re- 
sponse functions is an instance of the v-representability 
problem in density functional theories. 

In this paper, we introduce a stationary principle in 
TD DFT and use it to derive the KS equations for fre- 
quency dependent linear response. Our approach is to 
work in the time domain, and by turning on a harmonic 
perturbation Xv^ x ' (r) cosuit adiabatically, we induce a 
steady state linear response density An^^(r) coswi. We 
prove that the quasienergy (an analog of the Bloch 
quasimomentum for systems periodic in time), to sec- 
ond order in A, is a stationary functional of n^(r). 
If the linear response density Xn^ (r) cos wt can be re- 
produced by a KS system with an effective potential 
v s (r,t) = vi°\r) + Xv^^r) cos tut, then the station- 
ary principle implies that the exchange-correlation part 
of w^^r) is the functional derivative of the exchange- 
correlation part of the second-order quasienergy. We 
also find an analogous stationary principle and linear re- 
sponse KS scheme in the time dependent extension of 
one-matrix functional theory^ (1MFT). In 1MFT, the 
external potential can be nonlocal in space and spin co- 
ordinates, i.e., it acts as an integral operator with the 
kernel v(ra, rV). The corresponding many-body opera- 
tor is V = J2a*> I d 3 rd 3 r'^l{T)v{ra, rV)Vv(r'). Signif- 
icantly, a theorem analogous to the RG theorem has not 
been found in TD 1MFT. 

TD DFT for the special case of time-periodic exter- 
nal potentials has been studied previously . 15 ' 16 However, 
the scope of these and later approaches , 17 ' 18 which are 
based on a Hohenberg-Kohn-type minimum principle, arc 
severely limite d 19 ' 20 ' 21 because (i) the minimum principle 
is generally valid only for periodic potentials that have no 
Fourier component of frequency greater than the first ex- 
citation energy and (ii) one must assume the existence of 
Floquet states (reviewed below). The approach we pur- 
sue here is distinct because it employs adiabatic ramping 
in real time and relies on a stationary principle instead 
of a minimum principle; hence, our results are valid for 
all frequencies (except, of course, resonance frequencies), 
and we need not assume the existence of Floquet states. 

The paper is organized as follows. In Section [Til we 
review the basic properties of Floquet states, which are 
the fundamental states of a system with a time-periodic 
Hamiltonian. The Frenkel variation principle is used to 
derive a stationary principle for the quasienergy of Flo- 
quet states obtained from adiabatic ramping. In Section 
IIII1 we repeat the standard derivation of the frequency 
dependent linear response function, except that we em- 
ploy an arbitrary adiabatic ramping function. Theorems 
that establish a stationary principle and linear response 
KS scheme in TD DFT are proved in Section [TV] Anal- 
ogous theorems in TD 1MFT are proved in Section [V] 
A simple illustration of the KS scheme in 1MFT is pre- 
sented in SectionfVTJ We comment briefly on the question 
of the existence of / xc (w) in Section IVlIl 



II. STATIONARY PRINCIPLE FOR THE 
QUASIENERGY 

In this section, we derive a stationary principle for the 
quasienergy of Floquet states obtained by turning on adi- 
abatically a time-periodic perturbation of the external 
potential. 

We begin by reviewing the properties of the wave func- 
tion when the Hamiltonian is periodic in time, H(t+T) = 
H(t). If H{t) is an operator on a finite-dimensional 
Hilbert space, the Floquet theorem asserts that there ex- 
ists a complete set of solutions of the for m 22 ' 23 



* n (t)=(»(t)e- 



Ut + T) =£„(*), (5) 



which are called Floquet states or quasienergy eigen- 
states. The quasienergy, e„, is defined modulo 2ir/T. 
The periodic factor £ n (i) satisfies the equation 



{H(t)-id t )Ut) = e n (n(t). 



(6) 



If the system under consideration has an infinite- 
dimensional Hilbert space, the existence of Floquet states 
is not guaranteed by the Floquet theorem. For many sys- 
tems, it may be the case that H (t) — idt has no nontrivial 
eigenfunctions so that Floquet states do not exists 

Now consider an iV-electron system that starts in a 
ground state at t — — oo and experiences an adiabat- 
ically ramped (AR) periodic perturbation of the form 

Xv [ r ] {Y,t) = Xf{t/T)v^(r,t), where w«(r,i + T) = 
v^(r,t) and / = f(t/r) is a ramping function with 
time scale r. The many-body Hamiltonian is H T (t) = 
HW + XVP(t) with #(°) = f + W + V(°\ where f is 
the kinetic energy operator, W is the electron-electron in- 
teraction and = J d 3 r jj'^frj^frj^r). The ramp- 
ing function is an arbitrary smooth function that satis- 
fies /(— oo) = and /(oo) = 1. An example of a suit- 
able ramping function is (1 + tanh(</r))/2. The pre- 
cise form of / is inconsequential and will be left unspec- 
ified in the following. The functions v^\r, t) form a 
one-parameter family of perturbations, and the action 
of an ideal AR perturbation is realized by taking the 
limit t — ► oo at the end of the calculation. Although 
H T (t) is not exactly periodic, it is still possible to define 
a quasienergy if the system approaches a steady state in 
the limit (t, t) — » (oo, oo). 

Following Ref. [2a, we factor the wave function as 



*(t) =&-(*) exp i 



dt'K T (t')), (7) 



where 



K T (t) 



(£r\H T (t) - i8t\S T ) 

(Cr|er) 



(8) 



In accordance with the terminology of Ref. |25|, the over- 
all phase factor will be called the secular phase, and the 
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factor £ T = £ T (t) will be called the nonsecular wave func- 
tion. A system will be said to evolve adiabatically into a 
steady state if the following two conditions are satisfied: 
1) the nonsecular wave function tends to a unique func- 
tion £ = £(t) with period T in the limit (t, t) — * (oo, oo), 
i.e., if for all e > there exist t' and r' such that 
Ur{t) - < e when t > f and r > r' and 2) all of 

the electrons remain localized in a finite region of space 
for all time. Real systems, for which the spectrum usu- 
ally has a continuum component, are not expected to 
evolve into such a steady state due to the possibility of 
multiquantum ionization^ However, if the perturbation 
is harmonic and the driving frequency u> = 2ir/T is not 
a resonance frequency, the first-order term of the power 
series of £ T with respect to A evolves adiabatically into a 
unique harmonic function £W (see Section IlIII) . In such 
cases, the system will be said to evolve adiabatically into 
a first-order steady state. The periodic function £, if it 
exists, will be called the steady state nonsecular wave 
function, while will be called the first-order steady 
state nonsecular wave function. The quasienergy associ- 
ated with £ is, cf. ©, 



K = 



(e|H(t)-^^ t |c) 



(9) 



where H(t) is H T (t) without the ramping function /. 

If the system evolves adiabatically into a steady 
state, the Frenkel variation principle implies that the 
quasienergy is a stationary functional of For an ar- 
bitrary time dependent Hamiltonian H (t), the Frenkel 
variation principle states that W = is the solution 
of the Schrodinger equation with the initial condition 
* (to) = *i if 



(6W\H'{t) -id t \V) = 0, 



(10) 



where S 1 ^ = <5\l/(i) is an arbitrary variation that satisfies 
(t ) = 0. Setting H'(t) = H T {t) and adding to (ID]) 
its complex conjugate, we obtain 

8(V\H T (t) - idt\V) + idt{y\5V) = 0. (11) 

Substituting ([7]) gives 



(12) 



As we are interested in the response to an AR periodic 
perturbation, we now consider a one-parameter family 
of variations <5£ T (with parameter r) such that £ r + S£ T 
evolves adiabatically from the ground state to a steady 
state £ + <5£. Then, upon taking the following limit and 
time average^ the second term of (fT2")) vanishes, and we 
obtain 



lim — 

(t,T)->(oo,oo) T J t 



t+T 



which suggests the definition 

1 f t+T ,(^r\H T (t')-td t ,\^ T ) 

KJ& = lim - / dt' ^ 1 ; | , IW 

(t,T)^(0O,00) T J t \£r|£r) 



*+ T dt , {tj\H{t')-id t ,\C) 



(14) 



where the subscript v denotes the given «W(r,t), We 
define the domain of K v to be the space of all steady 
state nonsecular wave functions £ that can be obtained 
for some v^(r,t). It is possible to choose a larger do- 
main, but this is not necessary for our purposes. Let 
denote the steady state nonsecular wave function corre- 
sponding to the given v^'(r,t). A variation <5£ will be 
called an admissible variation if + 5£ is in the domain 
of K v . Thus, for every admissible variation there ex- 
ists a one-parameter family of variations <5£ T such that 
£ r + <5£ r evolves adiabatically into £ v + <5£. Hence, ([13")) 
implies the stationary principle 



SK v [£}=0 



(15) 



for an arbitrary admissible variation 5£ at £ = £„. In fact, 
this result follows from a direct calculation of the first 
variation, assuming only the existence of £„. We have 
carried out the above derivation based on the Frenkel 
variation principle because it will prove useful when Flo- 
quet states do not exist, a case to which we now turn. 

If the system does not evolve adiabatically into an ex- 
act steady state, the limit in (TT4"|) does not exist in gen- 
eral. Nevertheless, for an AR nonresonant harmonic per- 
turbation \v^{r)f(t/T) cosujt, the system evolves adi- 
abatically into a first-order steady state described by 
^q ' + A^ 1 - 1 , where ^ °^ is the unperturbed ground state. 
We now show that the second-order quasienergy is a sta- 
tionary functional of 

Consider the trial function 



(16) 



where r — (r-y, . . . , rjv) and £^(r, t) is an arbitrary har- 
monic function with frequency uj subject to the constraint 
(£o°^ = 0- Expanding ([T2"j) to second order in A and 
repeating the steps leading to 115]) . one obtains the sta- 
tionary principle 25 



(17) 



where 



+ <e (0) |^ (1) (i')|e (1) )+cc.). (18) 



dt' 5K T (t>) = 0, (13) at = e (i) 



Eq. (|17[) applies for an arbitrary admissible variation 



The trial wave function in (fT 
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which is specified only up to first order in A, is sufficient 
to obtain the quasienergy through third order ^21 



(19) 



This is analogous to a well-known fact concerning the 
variational estimate of a ground state energy: the error 
in the energy is second order in the error of the trial wave 
function. The first-order term of the trial function in (|16D 
can be expressed as 



which leads to 



(i)/ 



= k%[(P]+k21[£ ) ]; 



(20) 



kIM } ] = (£>\hv>-4°>±«\& ) ) 



+ ($ ) \VV\£>)+C.c 



(21) 



where V^> is the many-body operator corresponding to 
«W(r). 

In the above derivation, we made no assumptions 
about the higher order terms of the power series of £ T 
with respect to A. The series truncated at order N, 



(22) 



may also approach a unique periodic function. In most 
physical systems, there will exist an integer N' such that 
the truncated series will cease to approach a periodic 
function if N > N' , owing to multiquantum resonances 
Nu> = E k — Eq. The analysis of this section can be 
extended up to the highest order for which there are no 
multiquantum resonances. 



III. PERTURBATION THEORY FOR 
ADIABATICALLY RAMPED HARMONIC 
PERTURBATIONS 

In order to show that the nonsecular wave function £ r , 
to first order in A, evolves adiabatically into a unique 
harmonic function, we repeat the standard calculation 
of frequency dependent linear response. We employ an 
arbitrary adiabatic ramping function / = f(t/r) instead 
of the usual ramping function e 71 * with r\ — > 0. Without 
loss of generality, we take the perturbation to be of the 
form Av| 1} (r,i) = Au (1) (r)/(t/r) coswi with ui > 0. The 
calculation will clarify the sense in which the steady state 
linear response is independent of the precise form of the 
ramping function in the adiabatic limit r — ► oo. 

We consider an ./V-electron system with the Hamilto- 
nian H T (t) = + XV^\t). The system is assumed 
to start in the ground state at t = — oo. We expand the 
wave function as 



*(*)=X)°*(*) e " 



(23) 



where and Ek are the eigenstates and eigenenergies 
of H(°> . The initial condition is a k {— oo) = Sf-o- To first 
order (for k ^ 0), 



-iV, 



(i) 



All 



dt'f(t' /T)cos(cut')e iUJkot ' , (24) 



where V^, = ($fc|V'W|$o) andwj-o = E k -E Q . Treating 
first the e luJt component of cosojt, we obtain 



dt'f(t'/T)e l(ul+UJko)t ' 



- --V { ' 



2 k0 



f(t'M 



1 



i{io + Lo k0 ) 



J(uj+u k0 )t' 



1 tf^) 



1 



dt' 



i(ui + uj k0 ) 



from integration by parts. The last term vanishes in the 
limit (t, t) — > (oo,oo) if u> =/= —uJ k o, which can be shown 
as follows. Let 



I(r) 



dt 



dt 



J(ui+u k0 )t 



ds 



(25) 



where s = t/r. According to the Ricmann-Lcbesgue 
lemma, I(r) vanishes in the limit r — * oo due to the 
rapidly oscillating phase factor. The lemma requires only 
that the condition ds\df(s)/ds\ < oo is satisfied. 
This determines the degree to which the ramping func- 
tion is arbitrary: the steady state linear response will 
be independent of the precise form of the ramping func- 
tion provided only that the later satisfies /(— oo) = 0, 
/(oo) = 1, and the condition above. This conclusion is 
independent of the details of the system. The same analy- 
sis applies for the e~ lut component of cosuit if w ^ Wfco^ 
Combining the results, we find that the a k (t) approach 
the functions^ 



V {1) r 
2 

i sin(wi) 



cos(wi) 
1 



1 



1 



UJ - UJ k0 

1 



LU +U} k0 



UJ - Wfeo UJ + LOW 



(26) 



Thus, the first-order term of the nonsecular wave function 
approaches the harmonic function 

tP\t) = Y t $\t>- iUkot *k (27) 

The corresponding linear response density is 

» (11 (r,t) = ^2Rc(4 1) (<)e- 1 — *) (* |ft(r)|**) 



cos(ujt) y Vi 



EC 

fc#0 



1 



LV-LV k0 LO + UJ k0 



x ($ |n(r)|$ fc ), 



(28) 
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where n(r) = I ijj> (r)^>(r) and we have used the fact that 
the can be chosen real in the present case. The 
frequency dependent linear response function is readily 
identified as 



X(r, r';w) = ^ 



fc#0 



1 



1 



LU - LU k0 10 + Wfco 



x<* |n(r)|$ fc )<* fe |n(r')|$o>- (29) 



IV. KOHN-SHAM SCHEME IN TIME 
DEPENDENT DENSITY FUNCTIONAL THEORY 

We now turn to the main results of the paper. In this 
section, we prove two theorems that establish a station- 
ary principle and KS scheme in TD DFT for the special 
case of harmonic perturbations. 

In the following, we assume that the frequency depen- 
dent density response function x( w ) is invertible.— This 
is not always true. For example, x( w ) will not be invert- 
ible if it has any null eigenvalues apart from the trivial 
one corresponding to an arbitrary constant shift of the 
external potential. Mearns and Kohr*^ have given an 
explicit example in which x( w ) has nontrivial null eigen- 
values; however, these occur only for particular isolated 
frequencies lu p . The theorems below can be extended to 
such cases by requiring that (r) and n^(r) are or- 
thogonal to all null "directions." We also assume that 
lu > is not a resonance frequency, i.e., lu ^ (E k — Eq). 

Consider an electron system that starts in a nonde- 
generate ground state at t = — oo and experiences an 
AR harmonic perturbation Xv^ (r) f(t/r) cos tut. The 
system evolves adiabatically from the ground state with 
density n^(r) into a first-order steady state with den- 
sity n(°)(r) + \riu\r) cos tot. Since, by assumption, x( w ) 
is invertible, the linear response density (r) deter- 
mines the perturbation vjp (r) up to an arbitrary con- 
stant. The perturbation, in turn, determines the first- 
order steady state nonsecular wave function ^\ as seen 

(i) 



(i) 



(2) 



(r). We may now state the first theorem. 
Theorem 1. — The second-order quasienergy K), 
is a functional of nL (r), and it satisfies the stationary 

(2) 

condition 5Ky = for an arbitrary admissible variation 
5n2~\r) at n£'(r) = n$ v (r). 

Proof. - The second-order quasienergy is a func- 
= Uu(r) by composition of 
Therefore, the domain of K v [n^ ] is 



tional of n 



(i) 



and 

the LR VR space. A variation 8n\p (r) is admissi- 
ble if n£}v(r) + Sn^ir) is LR VR. In Section HH it 



was shown that K^^ 1 ^] is stationary for an arbitrary 
admissible variation S^ 1 ^ at = £y . As £W is 
a functional of on the LR VR space, there ex- 
ists an admissible S^ 1 ^ corresponding to every admis- 



sible Sn[}'. Therefore, Ktf' [n^ '] is stationary because 

8K {2) = (8Ki 2) /8& 1 ))(SSW/6n$P)Sn£ ) , suppressing co- 
ordinate integrations, and the first factor is zero. 

This stationary principle can be used to derive an ex- 
pression for the exchange-correlation potential of a lin- 
ear response KS system. The linear response KS sys- 
tem is a noninteracting system that experiences the po- 
tential v a (r, t) — vi°\r) + Xvi X \r) f (t / t) cos tut and re- 
produces the first-order steady state density m°)(r) + 

Ani, 1 "* (r) cos tot of the interacting system. Such a system 
will exist if the following two u-representability condi- 
tions are satisfied. 

Condition (la). The ground state density of the in- 
teracting system is noninteracting v-representable (VR- 
N). This means that there exists a system of noninteract- 
ing electrons with a potential (r) such that the ground 
state is nondegenerate and reproduces the ground state 
density of the interacting system. 

Condition (lb). The frequency dependent response 
function Xs(w) of the noninteracting system in condition 
(la) is invertible on the LR VR space of the interacting 
system. 

Theorem 2. — If an interacting system satisfies con- 
ditions (la) and (lb), then the exchange-correlation part 



( 2 ) r^Wi 



of Vs'(r) is given by 



8K (2 } 
8ntp (r) 



(30) 



where Kxc [nB^] is the exchange-correlation part of the 
second-order quasienergy. 

Proof. - In analogy with Ref. y, the exchange- 
correlation part of the second-order quasienergy is de- 
fined as 

l-t+T 



in (P7|) . Therefore, £y> is a functional of ni, (r) on the 
space of linear response densities that can be obtained 
for some vjj (r). This space will be referred to as the lin- 
ear response w-representable (LR VR) space. Let n^} v (r) 
be the linear response density corresponding to the given wnere 



3 ,n£\r)nX ) (T') 



(i) 



d 6 rd 



(31) 



s 



(2) k«l- 



and 



1 



T 



t+T 



dt'(^\f + V^-E^-id t ,\^) 



t+T 



dt\zP\f + vM-EW-id t ,\tP) 



and ^g 1 ** is the first-order steady state nonsecular wave 



function of the KS system. The last term of (|3T[) sub- 

(2) 

tracts the Hartree contribution, K H '. The above func- 
tional are universal in the sense that they do not depend 
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on v^\r). However, they do depend on the ground state 
density and the value of uj, though this dependence will 
not be indicated explicitly. With these definitions, one 
finds id 2) = Sf ] + Jd 3 r v£\r)n£\r) + + 
Hence, the stationary condition of theorem 1 gives 



(2) 



SS ° 2) +«W(r)+4 1 )(r)+ U «(r) 



= 0, 



(32) 



where 

} (r) = SK^/Sn&ir). The KS system also 
evolves into a first-order steady state, and the stationary 
condition for its second-order quasienergy K„ = S 1 ^ + 
/ (f 3 ri)| 1) (r)ni 1 '(r) will be identical to (02]) if 

«W(r)= ? ;W(r)+««(r) + 4 1 c )(r). (33) 

Since x s (w) is invertible, Vs(r) is uniquely defined. 

Thus, i>ic?(r) is given by (|3"U|). This completes the proof. 

The steady state linear response density of the in- 
teracting system can be computed from the expression 
nW(r,t) = Y,i =1 0* (O) (r, t) + c.c, where ^ (0) (r) 
are the ground state KS orbitals and (f>\ (r, t) are ob- 
tained from first-order perturbation theory for the single- 
particle Schrodinger equation 

idtM^t) = (-5V 2 + «W(r))^(r,i) 

+ A»W(r)/(t/T)cos(wt)^(r,t) (34) 

with the initial condition 0j(r, — oo) = ^'(r). We re- 
mark that theorem 2 can be extended to cases where 
the ground state density is not VR-N but rather EVR- 
N (noninteracting ensemble v-representable) by following 
an approach analogous to the one taken in the next sec- 
tion. 

Theorem 2 implies that the exchange-correlation kernel 
can be calculated as 

f xc (r, r'; ui) = — m x -^— . (35) 

The stationary principle for the second-order quasienergy 
does not entail a causality paradox because its basic vari- 
able n^(r) is time independent. It is straightforward to 
derive the following formal expression: 

1 

2 



drn«(r) W «(r). 



(36) 



V. KOHN-SHAM SCHEME IN TIME 
DEPENDENT ONE-MATRIX FUNCTIONAL 
THEORY 

In this section, we generalize the theorems of the pre- 
vious section to TD 1MFT. 



Before proceeding, it will be helpful to review some 
basic results from static 1MFT. The defining feature of 
1MFT is that it has the capacity to treat systems in 
which the external potential is nonlocal with respect to 
the space and spin coordinates. Accordingly, the neces- 
sary basic variable is the one-matrix (one-body reduced 
density matrix), which is defined as 

7(2, x') = N / dx 2 ■ ■ ■ dxNp(x, X2,... xn',x', X2, ■ ■ ■ xn), 



(37) 

where x = (r, a), f dx = J2 a I ( ^ r , an d p = 
J2n w ™ \^n) (^n\ is the TV-electron density matrix with 
ensemble weights w n such that ^ n w n = 1. By ex- 
tending the Hohenberg-Kohn theorem, 1 Gilbert proved^ 
that 1) the one-matrix uniquely determines the ground 
state wave function and 2) there is a universal energy 
functional E v [7] that attains its minimum at the ground 
state one-matrix. There is also a KS scheme in 1MFT. 
From the stationary condition for the energy functional, 
Gilbert derived the equation 



dx'v s (x,x')4>i(x') = ei<t>i(x), (38) 



where v s (x, x') = v(x, x') + SW/Sj(x', x) and W = W[y] 
is the universal electron-electron interaction functional. 
This equation can be interpreted as the single-particle 
Schrodinger equation for the orbitals of a noninteracting 
system (the 1MFT KS system). The potential v s (x,x') 
is a functional of the one-matrix. The ground state one- 
matrix of the interacting system can be obtained by solv- 
ing self-consistently (|3"5|) together with 



(39) 



7(2,2:') =^2fi<f>i(x)<f>*(x'), 



where fa are occupation numbers that satisfy fi — 
N and < < 1. Generally, fractional occupa- 
tion numbers are required to reproduce the one-matrix 
of the interacting system, not only the values and 
1 as in DFT. This scheme was originally described as 
paradoxica l 14 i 26 ' 27 because the stationary condition im- 
plies that essentially all of the £j collapse to a single 
level. Therefore, it appeared that the single-particle 
Schrodinger equation would not define unique orbitals. 
However, when the occupation numbers are shifted 
slightly from their ground state values, the KS equations 
have a self-consistent solution for a one-matrix that is 
close to the ground state one-matrix and for which the 
degeneracy is lifted, 2 - Thus, the correct ground state or- 
bitals can be obtained in the limit that the occupation 
numbers approach their ground state values. The ground 
state orbitals, which are called natural orbitals, are the 
eigenf unctions of the ground state one-matrix, and the 
corresponding eigenvalues are the occupation numbers ^ 
As the occupation numbers are fractional, it is useful to 
interpret the KS system as adopting an ensemble state. 

In the time dependent version of 1MFT, a general KS 
scheme has not been found. Such a scheme should have 
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the capacity to treat systems in which the time depen- 
dent external potential is nonlocal with respect to the 
space and spin coordinates. In TD DFT, the existence 
of the KS scheme is implied by the RG theorem. But 
in TD 1MFT, it is not known whether there is a theo- 
rem as general as the RG theorem, i.e., for as general a 
class of time dependence. The Bogoliubov-Born-Green- 
Kirkwood-Yvon hierarchy provides an equation of motion 
for the one-matrix. However, this equation contains the 
two- matrix (two-body reduced density matrix), and it is 
not known whether the two-matrix is a universal func- 
tional of the one-matrix when the external potential is 
nonlocal. As we are interested in frequency dependent 
linear response, we shall narrow our attention to the spe- 
cial case of AR harmonic perturbations. In this case, 
there is a KS scheme in TD 1MFT. Theorems 3 and 4 
below are generalizations of theorems 1 and 2. 

We shall need to refer to a different linear re- 
sponse function. The one-matrix response function in 
the time domain is defined as x( x ii x 'i> ti\ x 2 , x' 2 , t 2 ) = 
S r y(xi,x' 1 ,ti)/5v(x2,x' 2 ,t2). The frequency dependent 
one-matrix response function will be denoted = 
x( x l; x 'i'i x 2, x' 2 ; to)- Similarly, the one-matrix response 
function of the KS system will be denoted Xs( w )- in or_ 
der to emphasize the analogy between this section and 
the previous section, some of the notations will be dupli- 
cated. 

As the KS system can be interpreted as adopting an 
ensemble state, our first step will be to prove a stationary 
principle for an ensemble-weighted quasienergy. 

Consider an electron system that experiences an AR 
local or nonlocal perturbation Xv^ (x,x')f(t/T) costot. 
Suppose that the system starts at t — — oo in the ensem- 
ble state = J2 n w n |<E> n ) ($ n |, where $„ are orthogo- 
nal stationary states, no two of which are in resonance, 

i.e., E m — E n ^ u> for all m and n. Let 7i 1 l(x, x') de- 
note the steady state linear response one-matrix corre- 
sponding to the given Vw (x, x'). Also, let denote 
the second-order quasienergy that would be obtained if 
the system were to start in the pure stationary state 
<f>„ instead of the mixed state . It is convenient to 
introduce a notation in which Hermitian functions of 
[x,x') are expressed with respect to the complete basis 
of ground state natural orbitals 4>f\x). For example, 

vjj\ij) = J dxdx' <t>*^°\x)v£ (x,x')<j)j (x 1 ). We may 
now state the stationary principle. 

Theorem 3. — If the frequency dependent response 
function of the ensemble is invertible for a given to, then 
the ensemble- weighted second-order quasienergy 

K^^Y.™^ ( 4 0) 

n 

is a functional of = Ju\ij) and satisfies the sta- 

(2) 

tionary condition 5Ky = for an arbitrary admissible 
variation at 7^ = 7^. 

Proof. - Let be the first-order steady state 



nonsecular wave function that would be obtained if the 
system were to start in the stationary state $ n at t — 
—00. The ensemble- weighted second-order quasienergy, 
for fixed w n , is a functional of 7^ because each K„ 
is a functional of and each is a functional of 
7 i 1) . Each k£ ] can be shown to be a functional of 
by repeating the arguments leading to (f2Tj) for a sys- 
tem initially in the state <!> n , assuming, as we have, that 
<!>„ is not in resonance with any of the other stationary 
states. Each $p is a functional of because 7^ deter- 
mines (up to a constant), which, in turn, determines 
■ The first variation of [yffl] can be expressed 
as 5Ki 2) = J2 n w n (5Ki 2) /S^mk 1 ') /S^S^, where 
the coordinate integrations have been suppressed. A 
variation is admissible if 7^ + 5"/^ is in the 

LR VR space of the ensemble. The first variation van- 
ishes for an arbitrary admissible variation Sjffl because 
SK^' / 5£n = for all n, which follows from a straight- 
forward extension of the arguments in Sec. [TTJ 

Now consider an electron system that starts in a nonde- 
generate ground state at t — —00 and experiences an AR 
local or nonlocal perturbation Xvu(x, x')f{t/r) cosut. 
If the following two conditions are satisfied, then there is 
a linear response KS system in TD 1MFT. 

Condition (2a). The ground state one-matrix is 
noninteracting ensemble w-representable (EVR-N). This 
means that there exists a system of noninteracting elec- 
trons with a potential Vg (ij) such that the ground state, 

which may be an ensemble state , reproduces the one- 
matrix of the interacting system. 

Condition (2b). The frequency dependent one- 
matrix response function Xs(w) of the noninteracting sys- 
tem in condition (2a) is invertible on the space of all 7^ 
that (i) are LR VR in the interacting system and (ii) have 
no diagonal and degenerate components. 

The diagonal components of are simply the linear 
response occupation numbers — (ii), while 

the degenerate components are j^j\jk) + ~/u\kj) and 
—ijw (Jk) + i^\kj), where <fij and are any pair 
of occupationally degenerate natural orbitals, i.e., nat- 
ural orbitals for which /j -* = fj,°\ The diagonal 
and degenerate components correspond to null eigen- 
functions of Xs( w )j so they are not LR VR in the KS 
system. Therefore, the appropriate basic variable for 
the linear response KS system is t4 = 7o> (2:, x') = 
^ij7u\ij)<f>4 (x)<f>* (x'), where the prime indicates 
that the diagonal and degenerate components are ex- 
cluded from the sum. In effect, 7^ describes the orbital 
degrees of freedom but not the occupation numbers. Sim- 
ilarly, let denote the projection of the given perturba- 
tion Vu to the nondiagonal and nondegenerate subspace. 
Also, let 7L 1 ]u be the linear response corresponding to vi, ^ . 

Theorem 4- — If an interacting system satisfies con- 
ditions (2a) and (2b), then its first-order steady state 
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one-matrix j^(ij) + Xjw (ij) cos cot can be reproduced 
by a KS system with the potential v s (ij,t) — Vs(ij) + 
\v^P (ij)f(t/r) sinwt. The contribution to Vg(ij) from 
the electron-electron interaction is given by 



obtained for some VtP(ij). In order for 8K\f' to van- 
ish for an arbitrary admissible variation, the expression 
in brackets in (|4"3")l must vanish for all "directions" ex- 
cept the diagonal and degenerate directions. Eq. (|4"3"|) 
is identical to the stationary condition for the ensemble- 
weighted second-order quasienergy of a KS system^ with 
the potential v s (ij,t) 
if 



(2) 



(ij) 



SK, 



(2) 



(41) 



where K^l [7w ] i s the interaction part of the second- 
order quasienergy. 

Proof. — The proof is analogous to the proof of the- 
orem 2. The existence of v s (ij, t) follows from conditions 
(2a) and (2b). According to condition (2a), the ground 
state one-matrix can be reproduced by a KS system in the 

ensemble state pi 0) = J2 n V l**,n) The $ s ,„ 

can be taken to be Slater determinants of N natural or- 
bitals. To show that wJ*- 1 -* (ij) is the functional derivative 
of a universal interaction functional, we first define 

n 

i 

x^lt + tiW -Ek 0) -id t ,\4>P), 

where i and i)^ are one-body operators and £i*„ is the 
first-order steady state nonsecular wave function that 
would be obtained if the KS system were to start in the 
pure state 3> s ,n- The contribution to the second-order 
quasienergy from the electron-electron interaction is 



1 

T 



t+T 



dt'(^\w\^) 



(42) 



(2) 

We have not partitioned K- n l into Hartree and exchange- 
correlation terms because the linear response density 
riu \ which appears in the Hartree term, cannot be ex- 
pressed in terms of 7"£ alone, for it depends also on diag- 
onal and degenerate components of 7^. In terms of the 
above functionals, the second-order quasienergy can be 

written = S^ 2) + £^ v£ ] (ififfl (ji) + K$ t . Hence, 
the stationary condition (theorem 3) for the interacting 
system is 



6KP = 



(2) 



(43) 



for an arbitrary admissible variation S'fu (ij). A varia- 
tion Sjffl (ij) is admissible ii 7 fu; v (ij) + fi*Jw(ij) can be 



vf\ij) + Xv^(ij)f(tlT)smtJt, 



t(1) 



(ij) 



T (i) 



(ij) 



(ij)- 



(44) 



We remark that the KS perturbation in TD 1MFT 
must be advanced by a phase of ir/2 with respect 
to the given perturbation because the linear response 

\jw\ij) cosujt of the KS system has a phase delay. As 
the KS system can be interpreted as adopting an ensem- 
ble state, its one-matrix is governed by the equation of 
motion 



id t j(t) = [i+v.(t),y(t)] 

which gives, to first order in A, 



i(<t>P\i>?Ht)) = (<f>. 



(0) i-W/.ni ,(o)\ 

lv s (t)m ), 



for any pair of natural orbitals for which ^= /j . 
Thus, the steady state linear response of the interacting 
system can be computed from the expression 

iVut) = /i 0) <4 03 |# ) (*)>+/i 0) <4 1) (*)l4 0) ) 

;(4 0) -/i 0) )^(ii)cosa;t (47) 



(45) 



(46) 



p(0) 



1 



While the linear response KS scheme does not give the 
diagonal and degenerate components of 7^ (ij), they can 
be obtained instead by finding the stationary point of 
K^Ij^]^ This is analogous to the situation in static 
1MFT, where the occupation numbers, which are not 
determined directly by the KS equations, can be obtained 
from the minimization of the energy^ 

The linear response KS scheme implies the Dyson-like 
equation 



X(ij;kl;u) = x s (i?',W;u;) + x s (ij, 

mnpq 

x A(mn,pq;uj)x(pq,kl;uj), 



run: to 



(48) 



where A(w) = Sw^/Sju- Excitation energies can be 
calculated from the poles of the response function by the 
method proposed in Ref. [ill However, if the potential 
is local, it may be preferable to use rather than l|48p 
because the single- particle eigenvalues of the DFT KS 
system are often good approximations to the exact low- 
lying spectrum, while the 1MFT KS system provides no 
approximation at all due to its total degeneracy. If the 
potential is nonlocal, the inverse response function, which 
also contains information about the excitations, can be 
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(2) 

obtained from the second functional derivative of iQ , 
as seen from the following expression: 

K^ 2) = -- J dxidx' 1 dx 2 dx' 2 ^(x' 1 ,xi) 

X -1 (xi , x[ ; x 2 , x' 2 ; u)^ (x 2 , x' 2 ) 

+ f dxidx^vgXx^ati-yWfax!). (49) 



VI. ILLUSTRATIVE EXAMPLE 

In the previous section, it was found that the time de- 
pendent KS scheme in 1MFT does not directly determine 
the linear response of the occupation numbers. To clar- 
ify this aspect of the theory, we use the KS scheme to 
calculate the linear response in a simple example. 

We consider a simple version of the Hubbard model.— 
The electrons are confined to a discrete lattice, each site 
of which can accommodate up to two electrons. The 
electron-electron interaction is modeled by an on-site in- 
teraction. As a further simplification, we consider that 
there are only two sites and only two electrons^ The 
unperturbed Hamiltonian is 



£(0) 



-*E 



T C2o 



+ U (nif nil + ^2T«2|) + V^°\ 



(50) 



where c ia and Ci a are the creation and annihilation op- 
erators of an electron at site i with spin a and hi — 
J2a c l(7 C i(T- The first term of the Hamiltonian represents 
the kinetic energy by introducing "hopping" between the 
sites with energy parameter t. 

For = 0, the ground state is |$ ) = 

I +^4t c 2| +£<4t c U +2/<4t c 2J.) |0), where 
x = cos(7r/4 — a>o/2) and y — sin(7r/4 — ao/2) with 
tan ckq = U / At. 

Wc arc interested in the linear response to a non- 
local perturbation. For simplicity we shall consider 
only spin independent perturbations, so there will be 
only spatial nonlocality. Thus, the AR perturbation is 
XV (1) f(t/T) cosut with 



• (1) ct C 



= E4 w 

ija 

= E^~ 



(51) 



where we have introduced Pauli operators, e.g., a x = 
Sct( c 1o- c 2ct + c 2a .ci a ). The spatial one-matrix of a general 
state * is defined as 

7fe')-E( Vl 'lsW|*), (52) 

a 

which can be expressed in terms of the natural orbitals 
as 

7(tf) = E^*W#0')- (53) 



There are only two natural orbitals, 

cos((9/2)e-"^ 2 
sin(6»/2)e^/ 2 

sin(6»/2)e-^/ 2 
-cos(0/2)e^/ 2 



and 



(54) 



As a Hermitian 2x2 matrix, the spatial one-matrix can 
be expressed as 

7 = I + A (sin 8 cos tpa x + sin 8 sin ipa y + cos 9a z ) 
= I + j-a; 7 = (7 K ,7 y ,7z); (55) 

where A = (f a — fa) /2 and a is the vector of Pauli ma- 
trices. It is also convenient to express the one-matrix 
response function x( w ) with respect to the Pauli basis, 
e.g., Xxy — $lx/5v y . For the ground state $o ; we obtain 



xM 





V 



y 



2 >2 











OJ^ — U). 
^20 



5- y 



where LOko = E^ — Eq. The Kohn-Sham response function 





4A 



i 
-i 



has one null vector corresponding to the null linear re- 
sponse of the occupation numbers to a "diagonal" per- 
turbation Sv \cf> a ) (<f> a \ - Sv \4>b) (cf>b\. 

For ^ 0, the KS response function becomes 





-i cos( 



i cost 




i sin (9 sin -0 — i sin # cos?/' 



—is'm9 s'mtp 
i sin 8 cos ip 




This has the null vector (sin 8 cos ip, sin 8 sinip, cos 8), 
which is just the unit vector with polar angle 9 and az- 
imuthal angle ip- We observe that this vector is parallel 
to the 7 of the ground state, cf. (|55|) . which implies that 
any perturbation of the KS system, even a nonlocal per- 
turbation, can change only the direction of 7 and not 
its magnitude. The magnitude I7I = A is related to the 
difference of the occupation numbers (the sum is fixed, 
fa + fb = 2). As noted in Section [V] the occupation 
numbers of the KS orbitals are not changed, to first or- 
der, by any perturbation to the KS system. This is a 
general feature of the KS scheme in 1MFT. Neverthe- 
less, the linear response of the occupation numbers can 

(2) 

be obtained from the stationary condition SKi, = 0. 



VII. ON THE EXISTENCE OF f xc (cu) 

In this section, we explain why the invertibility of 
X(t, t'), which is established by the RG theorem, does 
not imply the invertibility of x( w ) f° r a pure frequency 
component. 
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The RG theorem implies that the inverse response 
function x" 1 ^^') is defined on the space Af(^ to ,to), 
which consists of all n^\r, t) that can be realized for 
the given initial state 'I' (to) = ^t by some perturbation 
that is analytic at t — to- In order to obtain x~ 1 ( UJ ) 
from the Fourier transform of X (t, t'); we m ust have 
X _1 (t,t') = x _1 (i - t 1 ). This will be the case only if 
the system is in the ground state (or a stationary state) 
of the unperturbed Hamiltonian at t = to- Therefore, 
the relevant space is J\f(^> gs , t ), where ^ gs is the ground 
state. The time t is arbitrary but finite. Thus, it can be 
shown that the invertibility of x(t— t') implies the invert- 
ibility of x(w) on the space J^f u (^gs,to), which consists 
of all nW(r,w) that are the Fourier transform of some 
nW(r, t) e Af(^> gs ,t a ). However, A/L(* ffs , t ) is too small 
to establish the invertibility of x( w ) f° r a pure frequency 
component to. In other words, it does not contain the 
elements nW(r)(<J(a; — f2) + 8(uj + Q))/2, corresponding 
to n^^(r,t) = (r) cos fit. Such elements are absent 
because a system in a perfect steady state with density 
n gs (r) + An^'(r) cos fit + 0(X 2 ) for all time is generally 
never in an instantaneous ground state. Therefore, there 
is no time to at which to specify the initial condition as 
required above. Hence, the invertibility of x( w ) and the 
existence of / IC (w) are not implied by the RG theorem. 

VIII. CONCLUSIONS 

One of the fundamental questions that can be asked 
about a quantum system is: How does it respond to a 
harmonic perturbation? The first-order response of the 
density to a weak perturbation is described by the lin- 
ear response function x(r,r';w). Considerable effort has 



been devoted to calculating x(r, r';oj) with the Dyson- 
like equation Q. This equation contains the exchange- 
correlation kernel f xc {u>). In this paper, we have shown 
that yr/ju) = 6 2 Ki 2 J/{6n^{r)6n^(r')), where 
K^xc [tIu ] is a universal functional. 

The RG theorem establishes the existence of time de- 
pendent KS equations, but it has not been possible to 
derive the exchange-correlation potential from a station- 
ary principle. The quantum mechanical action princi- 
ple does not provide a suitable stationary principle be- 
cause its density-functional formulation contains bound- 
ary terms<£ For the special case of harmonic pertur- 
bations, we have found a stationary principle for the 
quasienergy that can be used to derive the first-order 
exchange-correlation potential. 

If the external potential of a time dependent system 
is nonlocal, then it is not known whether a KS scheme 
exists in general. Although 1MFT has the scope to 
treat nonlocal potentials, a theorem as general as the 
RG theorem has not been found in TD 1MFT. By ex- 
tending the stationary principle for the quasienergy to 
TD 1MFT, we have shown that there is a KS scheme 
for the linear response of the natural orbitals to a har- 
monic perturbation. The KS system experiences an adi- 
abatically ramped perturbation of the form v s (x, x' ,t) = 
Vs(x,x') + Xvi [x, x') sinwt. The part of vi {x, x') 
due to interactions can be calculated from the functional 
derivative of a universal functional. In contrast to the 
DFT KS system, the linear response of the 1MFT KS 
system has a phase delay of 7r/2, so that the KS poten- 
tial must be advanced by a phase of 7r/2 with respect to 
the given external potential. 
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In Refs. @ and 0, generalized action functionals were de- 
fined, but they were not required to be stationary. 
In order to distinguish Fourier transform pairs such as 
x(r, r'; t — t') and \{r, r';ui), we will use the following no- 
tations: xH = X(r, r'-uj) and x(* - t') = x(r, r'; t-f). 
The double limit lim( t T )^( 00 00 ) must be taken with care. 
It can be evaluated as lim-r^oo limsup^^. 
This is a special case of the so-called 2n + 1 theorem. 
More precisely, we assume there exists an e > such 
that none of the eigenvalues Ek lie in the interval 
[E +u> — e,Eo + u + e]. 

The precise statement of the approach is 
lim r _ 00 limsup^^ limA^o |(afe(t) - Aa[. 1) (f))/A| = 0, 
where the order of the limits is important. 



As x( w ) i s an operator between infinite dimensional func- 
tion spaces, its invertibility can only be established with 
respect to these spaces. We do not pursue this issue here, 
and we shall simply assume that x is invertible on a suffi- 
ciently large space of linear response densities. 
The condition that the states of the ensemble are mutually 
nonresonant is automatically satisfied if /| ' 7^ 0, 1 for all 
i, which is the typical case, due to the total degeneracy of 
the KS system. 

If any occupation number is equal to or 1, then x( w ) 
will have a corresponding null vector. 28 Although a con- 
dition of theorem 3, namely that x( w ) i s invertible, would 
therefore be violated, the stationary principle remains valid 
for variations in the complement of the null space of x( w )- 
The static version of this model was treated by 1MFT in 
Ref . [H. 



